High-order Dy multipole motifs observed in DyB2C2 with resonant soft x-ray Bragg 

diffraction 
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Resonant soft x-ray Bragg diffraction at the Dy M4^5 edges has been exploited to study Dy 
muhipole motifs in DyB2C2. Our results are explained introducing the intra-atomic quadrupolar 
interaction between the core 3d and valence 4/ shell. This allows us to determine for the first time 
higher order multipole moments of dysprosium 4/ electrons and to draw their precise charge density. 
The Dy hexadecapole and hexacontatetrapole moment have been estimated at —20% and -)-30% of 
the quadrupolar moment, respectively. No evidence for the lock- in of the orbitals at Tjv has been 
observed, in contrast to earlier suggestions. The multipolar interaction and the structural transition 
cooperate along c but they compete in the basal plane explaining the canted structure along [110]. 

PACS numbers: 71.70.Ch; 75.40.Cx; 78.70.Ck 



Resonant x-ray Bragg scattering (RXS) has recently 
become a powerful tool in modern solid state physics 
to investigate magnetic, orbital and charge order phe- 
nomena associated with electronic degrees of freedom. 
Especially in transition metal oxides, and, possibly, ac- 
tinide compounds, various phenomena can occur simul- 
taneously and it is not easy to identify the exact role and 
interplay of the various order parameters. However, in 
rare earth based compounds, orbital and magnetic order 
occur more autonomously from the lattice as the 4/ elec- 
trons are shielded by the outer 5c? electrons. Electronic 
orbital motifs are usually labeled antiferroquadrupolar 
(AFQ) or ferroquadrupolar (FQ) but higher order multi- 
poles are not established by observations. Argon atoms 
in an excited state with ^=3 show a significant presence 
of the higher order multipole moments as concluded from 
electron scattering but in the solid state multipoles 
beyond rank 4 have never been detected. Yet the align- 
ment of higher multipoles to the surrounding charge in 
the lattice is likely. For 4/ electrons (/=3) multipoles 
up to rank 6 are present and this results in a undulated 
and aspheric charge density. Which multipoles dominate 
these orbital order transitions remains a rather contro- 
versial topic as exemplified in CeBg % NpOa 3] and 
URu2Si2 4]. Using soft x-ray resonant Bragg scatter- 
ing (SXRS), we have obtained the first direct evidence of 
high-order Dy multipole moment motifs in DyB2C2. 

DyB2C2 has attracted much attention lately as its high 
AFQ ordering temperature Tq=24.7 K allows to study 
this phenomenon conveniently. At room temperature, 
DyB2C2 crystallizes in the tetragonal PA/mbrn structure 
and undergoes a structural transition with small alter- 
nating shifts of pairs of B and C atoms along c at Tq 



(5] which reduces the symmetry to PA2/mnm 6]- Below 
T/v=15.3 K antiferromagnetic order (AFM) is observed 
with complex moment orientations due to the underlying 
orbital interaction J] . RXS at the Dy L3 edge has been 
used to elucidate Dy dipole and quadrupole motifs 0, 1^ . 
A dipole transition (El) occurs between Dy 2p and 5d 
shells and a quadrupole transition (E2) between Dy 2p 
and 4/ shells. The first and dominant process probes the 
quadrupolar order of the 5d states and the latter transi- 
tion probes the quadrupolar order of the 4/ states. The 
quadrupolar origin of the reflection has been confirmed 
by azimuthal scans (rotation around the Bragg wave vec- 
tor). 

Ground state wave functions of the Dy 4/ shell have 
been proposed based on the observed absence of E2 in- 
tensity at the (00 1 ) charge forbidden reflection which 
implied cancellation of quadrupole (rank 2) and hexade- 
capole (rank 4) contributions [9j. These wave-functions 
are consistent with a recent inelastic neutron scattering 
study where the orbital fluctuation timescale has been 
determined 1^ . An alternative explanation for the RXS 
data suggest that the hexadecapole moment is zero and 
the overlap between the El and E2 resonances cancels 
the E2 intensity However, this analysis is not sound 
[1^ . SXRS at the Dy M^^^ edges accesses the 4f shell 
directly with the El transition, without the complica- 
tion of overlapping El and E2 resonances. Our results 
show that the Coulomb (intra-atomic quadrupole) inter- 
action between the 3d and 4/ shells is significant. The 
ordered quadrupole moment of the 4/ shell in the in- 
termediate state reshapes the observed 3^5/2 and 3^3/2 
electron density and leads to a core hole splitting. It is 
shown that this phenomenon causes interference between 
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FIG. 1: Diffraction intensity of orbital ordering (OO5) reflec- 
tion in DyB2C2 taken with incident a and vr polarization at 
r=18 K (middle). Top graph shows the absorption recorded 
in electron yield and the bottom graph shows the diffracted 
intensity after correction for absorption. The inset shows the 
intensity as a function of temperature for i?=1282 eV and 
1291 eV. 

different pathways of the scattering amplitude allowing 
the extraction of high-order multipole moments of the 4/ 
shell that have rank 4 and rank 6 (hexacontatetrapole) . 

A DyB2C2 single crystal has been grown by Czochral- 
ski method using an arc-furnace with four electrodes and 
cut with (001) perpendicular to the sample surface. Sub- 
sequently, it was polished and aligned with 26° azimuthal 
angle. Zero degrees azimuth corresponds to alignment of 
the b axis in the scattering plane. The orbital order- 
ing (00^) reflection has been recorded at the Dy M4^5 
edges of DyB2C2 at the RESOXS end-station of the SIM 
beam-line at the Swiss Light Source. The Dy M4^5 ab- 
sorption edges have been characterized with fluorescence 
yield (FY) and electron yield (EY) at RESOXS and the 
BL25SU beam line at SPring-8, respectively. 

The energy profile of the charge forbidden (00^) reflec- 
tion, together with x-ray absorption (XAS) data (EY), 
is shown in Fig. ^ Multiple features are observed with 
much larger energy spread than the multiplet structure 
of the XAS data, in contrast with the single oscillator 
recorded at the Dy L3 edge jQj . A similar energy profile 
has been reported for the M5 edge of Ho at the magnetic 
(0,0, t) reflection but the central sharp feature (solid sym- 
bol in Fig.nj is absent The authors argue that the 
large photo-absorption causes this 'gap' at the edge 
because strong absorption is inherent to soft x-ray scat- 
tering. To characterize the width of the scattering 9/29- 
scans have been recorded for the (00 i) reflection at each 
energy. Comparison with FY and EY data shows part of 
the width is due to the limited penetration depth. Sub- 



sequently, the integrated intensity has been corrected for 
absorption as shown in the bottom panel of Fig. ^ Ab- 
sorption effects hugely reduce the intensity of the central 
feature but cannot account for the multiple features in 
the energy dependent intensity. The different character 
of those structures becomes even more apparent from the 
temperature dependence of the scattered intensity. The 
AFM transition is witnessed at 15 K as a gradual increase 
in intensity at 1282 eV, however this increase is absent 
at 1291 eV (see inset of Fig. This cannot be caused 
by absorption effects on a single oscillator resonance. It 
appears that the scattered intensity at 1291 eV is sen- 
sitive to the orbital order only and its gradual increase 
below T/v shows that there is no lock-in of the orbitals as 
pro posed from neutron diffraction and symmetry analysis 

The structure in the energy profile is similar for a and 
TT polarization of the incident x-rays and independent of 
the sample temperature in the AFQ phase. Azimuthal 
angles of 56° and 86° give the same spectral shape but 
show different overall intensity as expected from the az- 
imuthal angle dependence of a quadrupole. This illus- 
trates merely one atomic tensor is active in the scattering 
process. We propose that this particular shape of the en- 
ergy dependence of the (00^) reflection is caused by the 
splitting of the 3d core states which results in multiple 
interfering resonators and consequently adds structure to 
the energy profile. Previously a core hole splitting caused 
by the intra-atomic magnetic interaction between 5/ and 
3d was used to describe the unusual double Lorentzian 
energy profile at the M4 resonance in Np02 There 
is no ordered magnetism in the AFQ phase and we intro- 
duce the intra-atomic quadrupolar interaction to partially 
lift the core hole degeneracy. The energy splitting is de- 
termined by the strength of the intra-atomic interaction 
and we show that the relative amplitudes are determined 
by the 4/ wave function. Thus, the spectral shape is 
constant in the AFQ phase. Note that the diffracted in- 
tensity is caused by the bulk order and subject to critical 
fluctuations while the energy splitting is due to the lo- 
cal orbital moment of the 4/ shell and independent of 
its orientation. We will first encapsulate the resonant 
diffraction theory specific to our problem and introduce 
the intra-atomic quadrupole interaction to portray the 
observed energy dependence. 

The scattered photon intensity da/dO. is proportional 
to the squared modulus of the amplitude / [ig 

K,Q q 

which is described as a product of spherical tensors with 
rank K. with —K < q < K describes the electronic 
response of the sample. with —K < Q < K de- 

scribes the geometry of the experimental set-up and po- 
larization direction of the incident and reflected x-rays. 
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where R (r, x) are 
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FIG. 2: Observed energy profile (crosses) compared to 
fits of the theory including absorption correction with 
a) wave-function taken from Ref. (red), b) (tI) 
set to zero (dotted blue), and c) (Tf) as free pa- 
rameter (black). (rl) is normalized to 1, (Ti) —— 0.2, 
(r|)=-f-0.3, r5_^i=0.8 eV, r5_±3=2.7 eV, r5_±5=5.3 eV, 
Fa ,1=1.9 eV and Fa ,3=5.4 ev' 

2 ' 2 2 ' 2 

Lastly, DQ^{a, I3,j) rotates onto the coordinates of 
the experimental reference frame used for with Eu- 
ler angles a,P,"f. The Dy site symmetry of 2/m gives 
q=±2 and K—2 due to the absence of charge {K=0) and 
time-odd (magnetic) (1^=1) order at (00^) in the AFQ 
phase. 

describes the atomic resonant process which is 
commonly represented by a harmonic oscillator. In our 
particular case F^ is a sum of several oscillators at the 
El resonance created by the splitting of the core state. 
The amplitude {J, M) of each oscillator is labeled by 
the total angular momentum ^=f ,f and magnetic quan- 
tum number M of the core hole. 



rjA^{J,M) 



^ Aj-e(J,M)+ir 



(2) 



Here, E is the photon energy, Aj is the difference in en- 
ergy between the degenerate 3dj shell and the 4/ empty 
states and F j ^7/ the lifetime of the intermediate state. 
e( J, M)=[3M2 „ j( iy^Q_ jg ^Yve energy shift of the M 
core levels due to the intra-atomic quadrupole interaction 
Q J, similar as in Mossbauer spectroscopy. Qj is a prod- 
uct of the 3d quadrupole moment and the /-electron elec- 
tric field gradient experienced by the 3d electrons with 
Q±a/g±|=7/3. 

The amplitudes interfere and the branching ratio be- 
tween the two edges is defined as a real mixing parameter 
rj. {J, M) is constructed from the structure factor of 



the chemical unit cell '^^ 



<(J,M) = (-l)^"-^^^(2 
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r=0,l,...,2J, x=\K-r 
equal an even integer 



dictates that A2 



reduced matrix elements, 
..,|-ftr-|-r| and q+x and r+x 
The Dy site symmetry, 2/m, 
and 4':^=(r^) - (T^J where 
(Tg^) = (5fo|Tg^|5fn), is the atomic spherical tensor of 
the Dy 4/ shell and |5'o) its ground state wave function. 
It follows that A^{^,M) contains the hexadecapole 
moment (Tj) and hexacontatctrapolc moment (Tj) 
in addition to the quadrupole moment (Tj). In case 
e(J, M)=0 the terms proportional to (T^) and (T^) 
cancel and A^ (J) is proportional to (T|) as usual. 

Let us here summarize the consequences of this anal- 
ysis of the scattering. The overall intensity is well repre- 
sented by scattering from the quadrupole (T|) of the 4f 
shell and reflects the same azimuthal angle dependence 
as the quantity observed in the dipole transition at the 
L3 edge . The dependence on the higher 4/ multipoles 
arises from the splitting of the core hole states, which 
results in the M dependence of the amplitudes of the 
different harmonic oscillators (Eq. ^ . Correspondingly, 
these amplitudes are influenced by higher 4/ multipoles 
as shown in Eq. O (Wigner-Eckart theorem). Therefore, 
the core hole splitting allows, for the first time, to mea- 
sure the ordered hexacontatetrapole moment of an elec- 
tronic shell in solid state physics. Note, that this is not in 
contradiction with the fact that the integrated intensity 
of a dipole (El) transition is sensitive to tensors up to 
rank 2 (quadrupole), while a quadrupole (E2) transition 
is sensitive to tensors up to rank 4 (hexadecapole). 

Figure 12 shows the measured energy dependence com- 
pared to the above theory including absorption effect. 
This model describes the data well. Aj typically de- 
picts the onset of an absorption edge and ra is calcu- 
lated by Ref. 0] at 0.22 which is in excellent agree- 
ment with our SXRS data. Q5=— 0.4 eV and the res- 
onance positions are indicated in Fig. |2l The lincwidth 
increases progressively from 0.8 eV for M—±^ to 5.3 eV 
for M=±|. I ( |, ±|) I is relatively small. Variations 
in {T2) can be partially compensated by the line widths 
to yield a similar result. The best fit gives (Tj )=— 0.2 
and (T2 )=0.3 while the 4/ ground state wavefunction 
taken from Ref. 9] gives (r2')=-0.57 and (T|)=0.36 and 
describes the data also well. The corresponding fits as 
well as the effect of setting (T^ ) to zero are shown in 
Fig. [3 18] . The temperature dependence observed at 
1282 eV and 1291 eV supports our analysis since the res- 
onant cross section contains additional time odd terms 
due to the magnetic order below Tjy. These terms may 
add at e(|, ±|) while they cancel at e(|, ±i). 

The Dy 4/ charge density obtained from the best fit is 
illustrated in Fig. O It supports a 90° zig zag alignment 
along c as demonstrated by its protrusions at top and 
bottom. In the basal plane a 90° zig zag alignment along 
[110] is preferred so that the concave part points toward 
the convex part of its nearest neighbor. Yet, the absence 
of lock-in of the orbitals confirm that the orbital motif 
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FIG. 3: Dy charge density and muhipole motif in the AFQ 
phase of DyB2C2. The basal plane is indicated in gray. Note 
the 90° zig zag alignment of the Dy orbitals along c and the 
canted zig zag alignment along [110]. A spherical charge den- 
sity has been subtracted to emphasize the asphericity. 

in the AFQ phase is similar to that in the AFQ+AFM 
phase and the Dy charge densities are canted away from 
[110] JJ. Such an AFQ structure logically emerges when 
the pairwise movement of B and C ions is in competi- 
tion with the multipolar interaction and promotes paral- 
lel alignment of the orbitals along [110] [l9j. In contrast, 
the 90° zig zag alignment along c is supported by both 
mechanisms. This is consistent with a recent study of 
(Dy,Y)B2C2 where quasi one dimensional AFQ along the 
c axis has been reported [2(1 ]. 

We should add that our analysis in terms of core va- 
lence interaction is consistent with the multiplet struc- 
ture interpretation. In the M5 and M4 absorption the 
multiplet structure is ~10 eV and ~5 eV wide, respec- 
tively, which is mainly due to the strong 3d-4/ and 4/-4/ 
Coulomb and exchange interactions [l^]. The M^^^Af^^ 
and 3d^^2^f^^ final state multiplets contain 6006 and 
4004 levels with an average lifetime width of 0.27 eV and 
0.55 eV, respectively. The strong electron correlation 
leads to non-diagonal matrix elements in the quantum 
numbers LSJ of the core hole. The levels of the same M 
connected by these non-diagonal elements excite partly 
coherent. This makes it difficult to separate the core hole 
states. In the free ion we do not expect an energy shift 
for distributions of different \M\. This was confirmed 
for the x-ray absorption spectrum with Cowan's atomic 
Hartree-Fock code, where we found roughly an equally 
broad multiplet when switching off the quadrupole com- 
ponent of the orbital contribution of the core- valence in- 
teraction. This suggests that any energy shifts in the 
distributions will be due to multipolar ordering of the 4/ 
states. 

In this theoretical framework the core hole state is un- 
coupled from the 4/"+^ state and the effect of the 4/ 



multiplet structure is indirectly taken into account by 
the effective widths of the M distributions, which exceed 
the intrinsic life time width and describes the data sur- 
prisingly well. 

In conclusion, the Dy 4/ hexadecapole and hexacon- 
tatetrapole moments in DyB2C2 have been measured 
with SXRS and their magnitudes are determined at 
—20% and -f30% of the quadrupole moment, respec- 
tively. The orbital order remains unaffected by magnetic 
order below Tn and no lock-in of the orbitals takes place, 
in contrast to CeBg. The structural transition and mul- 
tipolar interactions cooperate along the c axis while they 
compete along [110]. These findings amply demonstrate 
a new extension of the resonant x-ray Bragg diffraction 
method for the observation of high-order electronic multi- 
pole motifs. Of particular interest are materials in which 
higher order multipoles are believed to be of importance, 
such as CeBg, URu2Si2 as well as the i?B2C2 and skut- 
terudite families. 
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